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On Oscillating Satellites. By H. C. Plummer, M.A. 

i. In his paper “ On Periodic Orbits” * (to be quoted here as 
Ch.) Professor Charlier has tried to obtain in a simple analytical 
manner some of the results already found by Professor Darwin in 
his well-known memoir f on this subject (to be quoted as D.). 
The orbits investigated by Professor Charlier are those of 
oscillating satellites— i.e. those small periodic orbits which can 
be described in the neighbourhood of the five points where exact 
solutions of the problem of three bodies are known to exist. 
Some of his results I have tried to discuss in a more general way 
in a paper J (to be quoted as PI.) in which I have considered the 
motion of a small body in the neighbourhood of centres of libra- 
tion, or points of zero force, in any rotating field of force. In 
the present paper it is sought to push the analysis beyond the 

* Meddclanden fran Lunds Astronomiska Obscrvatorium, No. 18. 
t Acta Mathcmatica , t. xxi. p. 99. 
t Monthly Notices, vol. Ixii. p. 6. 
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first degree of approximation previously obtained and to consider 
more particularly the orbits in the neighbourhood of the three 
collinear points of libration, especially those of families a and b r 
for which numerical results are to be found in D . 

2. As far as possible the notation of PI. will be retained. 
The equations of motion of the third mass in the most restricted 
form of the problem of three bodies, referred to axes rotating 
with the line joining the finite masses, are 

. 0a . 0a 

*-2y = ^,y + 2x = d - 

where 

’ 0- K''' + s) + .'( v+ s) 

The masses of the finite bodies are /a and v, and the unit of mass 
is such that y -\-v = i. The origin is taken at the centre of 
gravity of the two revolving bodies, these lying on the axis of x 
and their distance apart being the unit of length. The unit of 
time is such that the angular velocity of the axes is also unity. 
The partial differentials of a with respect to x and y are : of 
the first order, p and q ; of the second order, r, s , and t * ; of the 
third order, P, Q, R, and S. Then if the point (x, y) is a centre 
of libration, for which p = q = o, and the coordinates of the small 
mass are &+£, y+rj, the equations of motion may be written 

27) = r^+s// + i(P^ 2 4 - 2 Q^ + R/j 2 ) 4 - . . . 

V + 2Z 5 = s£+£> 7 +KQi 2 + 2 R £?7 + Sq 2 ) 4 - • • • 

so long as the body continues within the region where the expan¬ 
sions are legitimate. To the first order of the relative coordinates 
the solution is already known. It is 

S = h cos (pit —a), rj = k cos (mt—fi) 

where 

mA + (r + £—4) m 2 + rt—8 2 = o 
h 2 = a 2 (t + m 2 ), k 2 = a 2 (r+m 2 ) 
tan(a—/3) = — 2 m/s 

These equations, found in PZ., §§ 4 and 5+, leave the origin of time, 
depending on «, and the scale of the orbit, depending on a t 
arbitrary. 

3. To proceed to the second approximation;we may put 
£=/i,cos (ratf—a) + r, rj=k coz (rat— + 


* The symbol t also represents the time; but this can hardly cause any 
confusion. 

t The constant a used above is connected with the constant c used in PI. 
by the relation a 2 (4w 2 + s 2 ) = c 2 . 
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where h , k are small quantities of the first order, r/ of the 
second order. In virtue of the preceding relations among the 
constants, the equations of motion now become 

f-2r/=r" + ^ + U1 
?+ 2 ? = 8 k' + t n ' + V J 


(A) 


where if terms beyond the second order be neglected 

2lJ=:PA 2 cos 2 (mt— a) + 2QMcos (mt—a) cos ( mt—f 3 ) 

4 - R & 2 cos 2 (rat—ft) 

2V=Q li 2 cos 2 (mt—a) + 2 1 &hk cos (: mt—a) cos (mt — ft) 

+ S & 2 cos 2 (mt—ft) 

Hence 

JJ—a 1 -\-b 1 cos 2mt -\~c x sin 2 mt 
V—a 2 + b 2 cos 2 mt 4- c 2 sin 2 mt 

where 

4 cq=PA 2 4 2 Q/^ cos (a— / 3 ) + R & 2 
4&!=PA 2 cos 2a-f 2Q hk cos (a + / 3 ) + RA ; 2 cos 2/3 
4c I =P^ 2 sin 2a + 2Q^sin(a + / 3 ) + R& 2 sin 2/3 

and the corresponding equations for a 2y b 2y and c 2 are found by 

Then the equa- 


+ 61=0 
-f-Cj —o 
~\-b 2 —o 

qmht +sk 1 + (4m 2 + ^)^ 2 + c 2 =o 

From this latter set of four equations we obtain 

S. — (/^.rri 2 4 i)b x 4 * sb 2 4 ~ 4 ^c 2 

o.h 2 =. sb x —( 47 n 2 + r)b 2 — 4 i mc 1 

<).&!= —4m5 2 — (4m 2 4- t)c x 4- sc 2 

r.& 2 = 4 mb! +sc x —(4m 2 +r)c 2 

where 

t)=3(4ra4 4 -s 2 — rt) 

4. From these results obtained by proceeding to the second 
order of the dimensions of the relative orbit we have two 
general properties which characterise the orbits of oscillating 
satellites : 


substituting Q, R, S for P, Q, R respectively, 
tions (A) will be satisfied by the solution 

t! =f^4-h x cos 2 mt + k x sin 2mt 
r} f =f 2 -{-h 2 cos 2mt + k 2 sin 2 mt 

provided that 

*/i + S/a + ai=o 

Sfi + tf 2 + a 2~0 

and that 

( 4m 2 4- r)h j 4- sh 2 4- 4 mlc 2 

—4mh 2 4 - (4^ 2 4 - T)k x + sk 2 

sh 1 4- (4m 2 4- t)h 2 —4 mk x 
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(1) The orbits of any family are tautochronous to this order 


439 


of approximation. 

(2) The mean position of the small body is displaced by an 
amount of the second order from the corresponding centre of 
libration. 

In this way it should be possible to find, with little extra 
trouble, better approximations to the relative paths of oscillating 
satellites than the elliptic orbits discussed in Ch. The results 
given in D. for satellites of families a and b show that tauto- 
chronism ceases to hold at a comparatively small distance from a 
centre of libration. This fact restricts the application of (1) to 
very limited areas, and suggests that an independent study of the 
variation of the period in this class of orbit would be very 
valuable. 

5. The formulae are much simplified when the collinear centres 
of libration in the special problem are considered. In this case 
2/=o, and we have for the partial differentials of £1 





Q=o 

S=o 


The rule for the ambiguous signs is this : For the orbits between 


the finite masses (family a) we take -f/b — v ; for the orbits 
beyond v (family b ) we take + p, + v ; for the orbits beyond /u 
(family c ) we take —/*, — r . Now for comparison with D., 
P=tt> v== tt > an fi taking the values of r T and r 2 from D 
pp. 108-111, we obtain: 


Family a . Family b. Family c. 


r x 071751 1*34700 094693 

r 2 0*28249 0*34700 1*94693 

t 139 878 6*0955 3* 1660 

t - 5*4939 - i *5478 -0*0830 

P +65*076 —39* 2 79 +6*822 

R ~ 3 2 *538 +19*639 -3411 

m 26082 1*6763 10706 

> 6 8025 2*8099 1*1463 


I I 
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6. Further simplifications are effected by making a = 5 : 0 , /3==|-r. 
Then 

4« t = PA 2 + R& 2 = P(/i 2 —1& 2 ) = 4« 2 a I 
4^ = PA 2 —RA 2 = P(A 2 + ^A 2 ) = 4« 2 ^. n 2 
4c 2 = 2RAA = -r-Phk = 4« 2 ^. a 3 
1 ^2 . * 7^2 O 

and hence 

/i = —a*ajr 

hx — — (4m 2 +O a2 «3 +4mfl 2 a 3 
& 2 = 4ma 2 a 2 —(4m 2 -f r)a 2 « 3 

/ 2 = A 2 = Ic L = o 


The following numerical values result:— 

Family a. Family b. 

log 8 2-89534 2-09008 

Family c. 

1-21887 

log aj 

2'l6976„ 

1-49597 

0-27043, 

log Oa 

93843S 

9-65902* 

9-52077 

log a 3 

9-03335 

9 - 42739 » 

9-34366 

log/1/a 2 

I-0240I 

071096,, 

9-76992 

hxjd 2 

- 4*1352 

+ 2-6261 

-0-5487 

^ 2 /® 2 

-I*9209 

+ r580O 

-0-2895 

log hfa 

0-0584Q* 

0*05055, , 

001333, 

log hja 

O65893 

0-47482 

0-31736 


The coefficients h and A are taken of opposite sign in accord¬ 
ance with the equation (PL p. 9). 

(r+ra 2 )A sin (a-/ 3 ) = 2 mh 

7. Thus the equations of the three families of orbits are 
obtained. They are, for family a : 

l = +10*5685 a 2 — 1*1439 a cos 2*60822—4*1352 a 2 cos 5*21642 

*? — +4*5597 « sin 2*6082^ —1*9209 a 2 sin 5*2164^; 

for family b : 

£ = — 5'i4°° a 2 —1*1234 a cos 1*6763^ +2*6261 a 2 cos 3*35262 

= +2*9841 a sin 1*67632+1*5800 a 2 sin 3*35262; 

and for family c : 

£ = + °' 5&&7 ^ 2 —1*0312 a cos 1*07062—0*5487 a 2 cos 2*14122 
= +2*0766 a sin 1*07062—0*2895 a 2 sin 2*14122 

8. The next step is to connect the parameter a with the con- 
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stant C of the Jacobian integral Y 2 = 2H—C. The equations 
representing the orbits being of the form 

£ = h' cos mt . d -f (fi 4-A/60s 2 mt ). a 2 
r) ■= k r sin . a 4* A 2 ' sin 2mt . a 2 

we have at the time t = o, 

£o'= A r . fls + ( t /i / + A/) a- \ r \0 = 0 
£ 0 = o ; iy 0 = wtA'. «+ 2mA 2 ' . a 2 

Hence to the third order in a, 

Y 2 = m 2 k ' 2 . a 2 + 4m 2 A'A 2 '. a 3 

2^ = C 0 +r£ 0 2 -f-i Pf Q 3 

= C 0 +rA' 2 . a 2 + {2rh f (/ I f + A/) +£ PA' 3 } . «3 

where C Q is the value of 212 at the point of libration. There¬ 
fore 

C—C 0 = (rh ' 2 — m 2 k f 2 )a 2 -f- {2rA'(/ I f -f- hj )—4 m 2 k'k 2 '+i PA' 3 } «3 

Now the coefficient of a 3 herein vanishes, as can be proved 
without difficulty. A useful check on the accuracy of the 
numbers occurring in § 7 is thus obtained, and the last equation 
can therefore be written in the simple form 

C 0 — C = a* 

which is equivalent to a formula in PL (§ 7). The numbers for 
the three cases are these : 

(a) (b) : ’ (c) 

C 0 3*65292 3 * 534 i 8 3*17322 

m*—rt 123*13 i)*33 i*576 

9. That the coefficient of a 3 in the expansion of C—C Q must 
vanish is clear from the fact that its sign would be changed by 
taking the initial conditions corresponding to mt=zw. And the 
verification is easily obtained from the general relations of § 3 as 
simplified in §§ 5 and 6. For we have 

hf —zm t + m 2 , x 

k r + m 2 —2 m v ' 

4 Vfx + PA^ 2 4 - HA' 2 = o ... . (2) 

(4m 2 +r) A/+4mA 2 '+i(P/6 / 2 -RA' 2 ) = o ... (3) 

4mhi+(4m 2 + 1 ) A 2 '-f£ RA'A' = o ... ... (4) 

Let (3) be multiplied by m and (4) by (m 2 *f-r). The difference, 
when account is taken of (1), becomes 

2,mrh/ 4- 3 (m 2 +r) m 2 A 2 '—(PA' 2 + 3PA' 2 ) = o 

1 1 2 
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which again becomes when multiplied by 2 h'/$m 
2rh'h 1 '-4m*k% , -%h' (P7</ 2 + 3 RA;' 2 ) = o 
Hence by (2) 

2rh r (fi + hi) — 4m 2 k'k 2 '+±Ph '3 = o 
which is the result required. 

10. Four examples of orbits of family a are given in D . 
(pp. 233-238), the values of 11 C being 40 o, 39*0, 38*5, and 38*0. 
These give the numbers 

11 (C 0 —C) = 0*1821, 1*1821, 1*6821, 2*1821 
log a = 8*06428, 8*47046, 8*54705, 8*60356 

Whence the following four sets of equations for the orbits 
result : 

(a) £ = + *0014 — *01^3 cos mt —'0006 cos 2 mt 

t) = +*0529 sin mt —*0003 sin 2 mt 

(/ 3 ) £ = -j- *0092 — *0338 cos mt — *0036 cos 2 mt 

t) = -f*i347 sin mt —*0017 sin 2mt 

(y) £ = -f *0131 — *0403 cos mt —*0051 cos 2 mt 

7 ? = +*1607 sin mt— *0024 sin 2mt 

(l) £ = -f *0170 — *0459 cos mt —'0067 cos 2m£ 

77 = +*1830 sin *0031 sin 2mt 

For comparison with Professor Darwin’s calculations the fol¬ 
lowing numbers may suffice. They give the points where the 
orbits cross the axis of x and the coordinates of the points for 
which mt = |-rr; for the latter the ordinates corresponding to 
the computed abscissae are found by interpolation among the 
coordinates given in D . 


mt = o 

7 T 




f. 



^1 7 

v 1 

(a)--0125 

+ *0141 

+ -0016 

+ *0020 

+ ^529 

P.--OI25 

+ 0147 

+ *0022 

ft 

+ *0530 

($) — ‘0282 

+ *0394 

+ -0112 

+ 0128 

+ '1347 

£>.-•0304 

+ •0413 

+ *0109 

9 9 

+ *I 395 

(y) — *0323 

4 *0483 

+ *0160 

+ -0182 

T *1607 

£.-0361 

+ •0523 

+ *0162 

99 

+ *1691 

( 3 ) --0356 

+ 'O562 

4- 0206 

+ *0237 

+ •1830 

£>.-•041 

+ *c6o 

+ *019 


+ •197 


In every case £' is slightly less than the abscissa corresponding 
to the maximum ordinate in D. } as should be the case. 
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11. Two examples of orbits of family b are given in D. (pp. 239- 
242), the values of 11 C being 38*5 and 38*0. These give 

11 (C 0 —C) = 0*3760, 0*8760 

Jog a = 8*64750, 8*83116 

and the two corresponding pairs of equations for the orbits are 

(u) £ = — *oioi— *0499 cos mt + *0052 cos 2 mt 

t) — +*1325 sin mt -f *0031 sin 2 mt 

( f 3 ) £ = — *0236 —*0762 cos mt -f- *0121 cos 2 mt 

rj 555 + *2023 sin mt 4- *0073 sin 2 mt 

The results of a comparison with D . are given in a form 
similar to that used in the preceding paragraph. 


mf= o 

IT 



\ir 

f. 

(a)--0548 

ft 

+ 0450 

fo + ft 
— *0098 

-•0153 

V 

+ •1325 

■D.-- 055 ' 

+ *0445 

-•0106 

»» 

+ •1312 

N 

OO 

O 

ai 

+ •0647 

- -0230 

-•0357 

+ *2023 

D. — *0875 

+ *0679 

— •0196 

»» 

+ -1987 


In both cases £' is slightly greater than the abscissa corre¬ 
sponding to the maximum ordinate in Z>., as the above equations 
require. 

12. Oscillating satellites belonging to the family c have not 
been investigated by Professor Darwin, and there seems to be 
no means of judging the accuracy of the above results for this 
family. The comparisons contained in §§ to and 11 for satellites 
of the families a and b seem to show a fairly good agreement, 
especially, as was to be expected, in the case of the smaller 
orbits. For the larger orbits, particularly those of family a, the 
results are by no means so good ; and it may be noted that 
Professor Darwin himself found the orbit (b) “exceedingly 
troublesome.” But by going to a second approximation it is clear 
that a better notion has been obtained both of the mean position 
and shape of the orbit than was possible with the elliptic paths 
investigated by Professor Charlier. The results found in Ch, 
(e.g. p. 25) are really better than could fairly have been 
anticipated owing to the fact, mentioned in § 8, that the 
coefficient of a? in the development of C vanishes. Further 
progress in the study of the particular problem considered in the 
present paper seems to depend, as noticed in § 4, on a knowledge 
of the variation of the period and the effect of that variation on 
the form of the orbit. 

University Observatory , Oxford: 

1903 June 11. 
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